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ABSTRACT. The problems associated with measurements of star formation rates 
in galaxies are briefly reviewed, and specific models are presented for 
determinations of current star formation rates from Ha and FIR luminosities. 
The models are applied to a sample of optically blue irregular galaxies, and 
the results are discussed in terras of star forming histories. It appears 
likely that typical irregular galaxies are forming stars at nearly constant 
rates, although a few examples of systems with enhanced star forming activity 
are found among HII regions and luminous irregular galaxies. 


1 . INTRODUCTION 

In this paper we briefly review methods for deriving star formation rates 
(SFRs) in galaxies with an emphasis on techniques based on the Ha and IRAS 
infrared luminosities (see also Gusten and Mezger 1982). Conversion of 
measured quantities into SFRs involves model— dependent assumptions, and it is 
interesting to intercompare SFRs derived by different methods as a means to 
test the reliability of the models. The blue irregular galaxies which we 
discuss include objects that are espcially well suited for this test for a 
variety of intrinsic and practical reasons outlined in Hunter and Gallagher 
(1986). B 


The strong emission lines and blue colors of galaxies in our sample 
sometimes are taken as evidence for "star burst" events. This view conflicts 
with results from the analysis of Gallagher, Hunter and Tutukov (1984; 
hereafter GHT), where a near constant SFR was found to fit most blue Irrs in 
agreement with earlier studies by Searle, Sargent and Bagnuolo (1973). The 
IRAS data base provides a new way to explore the star burst issue (as 
emphasized in many talks and poster papers during this conference), and we 
have developed a working definition for star bursts that is applicable to most 
types of galaxies. 


2. MEASURING CURRENT STAR FORMATION RATES IN GALAXIES 

To illustrate the factors involved in calculating SFRs, consider some 
quantity Q which is a measurable property of a stellar population in a galaxy 
(i.e. Q could be a luminosity, color, or mass estimate). Formally we then 
have 


Q = I t !yf m 9(®) <Km) w(r) g(t) dm dV dt • H(D). (1) 

* See erratum, page 257. 
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This equation illustrates the many parameters that affect any apparently 
simple measurement of Q: q(m) is the weighting factor in stellar mass, 4>(ra) is 
the stellar initial mass function, g(t) is a description of the time evolution 
of the system, w(r) takes internal structural properties such as extinction 
into account, and~H(D) is the distance factor which is (4irD) 2 for fluxes. 

Thus we require either good models or simple galaxies if we are to obtain 
reliable measurements of Qs that can be converted into SFRs, and in our work 
we have chosen the latter route in concentrating on studies of irregular 
galaxies (see Gallagher and Hunter 1984; Hunter and Gallagher 1986). 


To get a current SFR for a galaxy, we must observe a Q that is sensitive 
to short-lived stars. If x is the time scale for significant variations in 
SFRs over galactic scales, then we require that q(m*) ■+■ 0, where the stellar 
evolution lifetime t for stars of mass m* is x(m*) ~ tg. Unfortunately, we 
do not currently have good theoretical predictions for x<j, and furthermore Tg 
may vary between galaxies. A reasonable choice is x^ < 10 yr, which g 

corresponds to typical dynamical time scales within galaxies. For Xg — 10 
yr, the critical stellar mass is m* ~ 5 M Q . On the other hand, one should not 
choose too narrow a time interval to probe SFRs, since statistical noise 
associated with small samples of objects or small numbers of star forming 
sites will become a problem. For example, counting the luminosities of HII 
regions in the short-lived compact phase (Habing and Israel 1979) would not 
necessarily give a statistically reliable indicator of the current SFR in a 
galaxy. 

To measure current galactic star formation rates, we therefore must 
estimate the numbers of OB stars, which is a well known result. A variety of 
measurable quantities Q can be used to trace OB stellar populations in 
galaxies: (1) Rocket ultraviolet luminosities reflect the numbers of hot stars 
in galaxies, but the interpretation is complicated by major extinction effects 
which are not easily disentangled (e.g. Huchra et al. 1983; Donas and 
Deharveng 1984; Lamb, Hunter and Gallagher 1986). (2) In a few nearby 

galaxies, OB stars can be directly counted and the resulting luminosity 
functions converted to star formation rates (Dennefeld and Tammann 1980; 
Hoessel and Danielson 1983). A difficulty in this approach lies in obtaining 
complete samples of OB stars, which is proving to be a major challenge even in 
the Magellanic Clouds. (3) Thermal radio emission and hydrogen recombination 
line luminosities measure the fluxes of Lyman continuum photons and thus the 
populations of hot, young stars. This approach has been successfully adopted 
by several groups including Smith, Biermann and Mezger (1978), GHT, Hunter and 
Gallagher (1986), Israel (1980), and Kennicutt (1983). These results also are 
sensitive to effects of dust and to the choice of model to convert observables 
into Lyman continuum luminosities (see the excellent summary by Glisten and 
Mezger 1982). (4) The far infrared luminosity may arise from dust heated 

primarily by OB stars. Our lack of understanding of the details of radiation 
processes responsible fcr L(FIR), as discussed by several papers during this 
meeting, introduce uncertainties into the choices of w(r) and q(m) that are 
needed to drive SFRs. In the pre -IRAS era, L(FIR) primarily was used to probe 
SFRs in 'star burst' galaxy candidates (e.g. Telesco and Harper 1980; Rieke et 
al. 1980; Gehrz, Sramek and Weedman 1983) and only recently has this Q been 
applied to galaxy populations in the general sense. (5) Radio continuum 
luminosities of galaxies are well-correlated with all other indicators of 
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current SFRs (see Klein 1982; Kennicutt 1983; de Jong et al. 1985; Helou et 
al. 1985; Sanders and Mirabell 1985), but suffer from a lack of understanding 
of the basic physical process which produces the correlation (i.e. both 
w(r) and q(m) are not properly known). Thus at present radio continuum data 
must be taken as a secondary method for deriving SFRs, which relies on other 
techniques for absolute calibration. 


3. A COMPARISON OF SFRS BASED ON FIR AND Ha LUMINOSITIES 

In the spirit of eq. (1), we compress the problem of measuring a 
corrected Ha luminosity L°(Ha) from an Ha flux corrected for Galactic 
extinction F°(Ha), 

L°(Ha) = F°(Ha) • exp x(Ha) • (1+y*) . (4xD 2 ). (2) 

Here D is the distance, y* is a correction factor for underlying Ha stellar 
absorption, and x(Ha) is some mean optical depth in the Ha line due to 
internal dust in the galaxy. From simple recombination theory, the Lyman 
continuum photon luminosity N c is proportional to L°(Ha) , N c = 7. 3x10* 1 L 0 ( Ha ) 


i ! i i ?^ ed „ t0 t l le stellar population through the initial mass function 
(IMF), and following the notation of GHT with a Salpeter IMF, 


N c -« c 


f m u , 

J„. ( " 


-2.35 


) [n (m)t(m)]ri dm. 
c 


(3) 


We have adopted the Renzini-Tinsley approach in assuming a constant SFR over 
the time interval of interest here, in which case the time evolution of the 
luminosity for stars of mass m, n c ( m ), is replaced by an average over the 
appropriate stellar lifetime at mass m for the production of H-ionizing 
photons (see Renzini 1981). The other parameter n is an efficiency factor, 
and we wish to solve for the normalized SFR, a (M(Ha) = 5.8a M yr with a 
minimum stellar mass of 0.1 M Q ; see GHT for details.) C ° 


The integral in eg. (3) is also sensitive to the choice of the maximum 

stellar mass m u , as n c (m)t(m) rises nearly as steeply as the Salpeter IMF 

falls to masses of ^100M (Glisten and Mezger 1982). For an m, = 100m-, we 
estimate u ° 


M(Ha) = (2.5xl0“ 8 ) n” 1 L° (Ha) M q yr -1 (4) 

where L°(Ha ) is measured in units of the sun's bolometric luminosity. The 
coefficient in eq. (4) depends on several factors, including the choice of IMF 
parameters, stellar evolution models, and the derivation of n ( m ) as a 
function of time, effective temperature, surface gravity, abundance, etc. 
Hidden in the efficiency factor n are the escape of ionizing photons and the 
loss of ionizing^ photons due to the effects of dust. Since most of the 
contribution to ft(Ha) must come from stars with lifetimes x(m*) < 10 7 yr; Ha 
observations in principle provide particularly clean snapshots of star forming 
processes in galaxies. 
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An SFR can be derived from measurements of L(IR) under the assumption 
that massive stars play a major role in heating dust. The analog to eq. (3) 
is then (cf. Hunter et al. 1986a) 

L*(FIR) = a 1 /™“ (m~ 2 * 35 ) ? (L(m) tCm))^ dm. (5) 

This is a more extreme simplification than in the Hci case, as the efficiency 
factor 6 should be a function of stellar temperature and the geometry of the 
galaxy should be taken into account in calculating the absorption of 
radiation. We will assume that 6 ~ 0.5 for all stars with an effective 
temperature above some threshold. The sum is over stellar evolutionary 
phases, and following the Renzini-Tinsley approach, we need only bolometric 
energies which stellar evolution models predict reasonably well (cf. Renzini 
1981). Finally, we must account for infrared luminosity that is not 
associated with young stars, e.g. we probably should exclude the cold dust 
luminosity component (see Hunter et al* 1986a). Thus for an observed L (IR), 
we have L*(FIR) = 6L°(IR) where 6 < 1. 

For m* = 10 M q and modern stellar evolution models which include mass 
loss (e.g. from Maeder 1981 ; Brunish and Truran 1982), we find 

Mj (FIR) = 2.5xlO“ 10 g -1 S L°(IR) M^r" 1 (6) 

where L°(IR) is again in bolometric solar units. The relevant time scale for 
eg. (5) is i(m*) ~ 2xl0 4 * * 7 * * yr, and we are including only stars and evolutionary 
phases with log T > 4.3. If cooler stars are effectively heating dust, then 
the coefficient in eq. (6) is reduced. For example, reduction of the stellar 
temperature limit to log T g > 4.0 corresponds to m* ~ 2 M Q and_|(m*2^~ 10 yr, 
while the coefficient changes only slightly to 2.1x10 M Q yr L Q * 

Taken on their own, the FIR data thus involve a time average of the SFR 
over a somewhat indeterminate period, which could vary depending on physical 
conditions within galaxies. Dusty galaxies can have high optical depths for 
the absorption of radiation at visible wavelengths, and even comparatively 
cool stars can contribute to L°(IR). Most blue irregular galaxies lie at the 
opposite extreme, and have low dust optical depths in the visible (Hunter and 
Gallagher 1986). In these types of galaxies, UV heatings by younger stars is 
likely to be favored, and thus L°(IR) may be a good indicator of the current 
SFR. 


4. APPLICATIONS TO REAL GALAXIES 

4.1 Current SFRs 

Blue irregular galaxies provide a way to test our ideas about SFRs and 

evolutionary processes in galaxies for the following reasons: (1) We have 
collected a large body of optical data, which includes Ha fluxes (cf. Hunter 

and Gallagher 1986). (2) Due to the low internal optical depths found in most 

Irrs, star forming sites usually are optically visible (Hunter 1982). (3) 

Irrs have the bluest optical colors for 'normal' (i.e. non-active) galaxies, 
and thus are logical candidates for star formation bursts on galactic 


170 



STAR FORMATION RATES IN BLUE GALAXIES 


scales* (4) The majority of the less luminous Irrs are structurally simple 
and close to single zone systems, which allows primitive models to be used in 
interpreting observations. 

We have recently collected IRAS observations of Irr galaxies included in 
our various optical samples (see Hunter and Gallagher 1986; Hunter et al* 
1986a, b). These data are shown in Figure 1, where the Irrs basically form a 
continuation of the normal spirals to higher far infrared color temperatures, 
and perhaps somewhat more surprisingly, to higher values of L(IR)/L(B). A 
complete discussion of these data is currently in preparation by Hunter Rice 
Gallagher and Gillett. 



log S(100)/S(60) 


Figure 1. Irr galaxies are plotted on a luminosity-color diagram after de 
Jong et. al. (1984). L(B) is the blue luminosity, and S(60) and S(100) are 
IR AS flux measurements at 60 and 100 microns. L(IR) is derived from vf(v) at 
80 microns. The various structural class samples of Irr galaxies defined by 
Hunter and Gallagher (1986) are shown, and the approximate area occupied by 
spirals is outlined. 


The two indicators for current SFRs, ft(Ha) and ft (IR), can be directly 
intercompared using the data for Irrs. Making the approximation that (1+y*) 
= 1.0 in eq. (2), we have 


M(Hct) 

M^IR) 


lf>A exp x(Ha)8L(Ha) 

100 tTS~ lTTrT * 


(7) 
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For the assumptions of high efficiencies, (6 = n - D and low T(tta), we obtain 
the limit that L(IR) = 100 L(Ha) if both luminosities measure the current SFR, 
and star formation has continued at a constant rate over > 2xl0 7 yr (cf. eqns. 
5 and 6). A minimum ratio of L(IR) = 16 L(Ha) can also be derived by assuming 
that, as in HII regions, equal amounts of power from 0 stars goes into dust 
heating and the eventual production of Ha emission (see Hunter et al. 

1986a). The minimum model thus is appropriate where L(IR) is dominated by 
contributions from HII regions. The results of these models are shown in 
Figure 2. 



Figure 2. Ha and FIR luminosities are plotted for Irr galaxies (same symbols 
as in Figure 1) and spiral galaxies for which Ha fluxes have been measured by 
Kennicutt and Kent (1983) (plotted as stars). Model relationships are also 
shown. 

Figure 2 is encouraging in that L(IR) ot L(Ha) over a factor of ~10^ in 
luminosity and the proportionality constant lies near the expected value of 
100. Thus we see tht L(IR) empirically is a useful tracer of the current 
SFR. This diagram also shows that many details remain to be sorted out. 
Luminous spirals and Irrs have high L(IR)/L(Ha) ratios, which could be due to 
high optical depths or heating by older stellar populations. Just as puzzling 
are the normal Irrs where L(IR)/L(Hct) falls near our minimum estimate. It is 
possible that only very massive stars are contributing to the observed L(IR) 
in these systems, and further work is needed. 
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4.2 Star Formation Histories and the Existence of Star Formation Bursts 

Following the approach of GHT, we can chart the evolutionary histories of 
galaxies by comparing SFR indicators that are sensitive to long term integrals 
over the SFR with estimates of current SFRs. GHT and Kennicutt (1983) have 
argued that most disk galaxies with active star formation have produced new 
stars at roughly time-constant rates during the past several Gyr. We can test 
the constant SFR, fixed IMF model by using the blue luminosity L(B) to 
estimate the average SFR over ~3Gyr (GHT). 

In this case, a simple model which includes evolved stars as in eq. (5) 
and adopts a Salpeter IMF predicts L(B)/L(IR) ~1, where we have taken a mean 
optical depth at B of Tj » 1.0. Note that L(B) is the luminosity in terms of 
solar luminosities on the Johnson B system (one must be careful to properly 
define luminosities in the various color bands!). The results of this model 
are shown in Figure 3. Most galactic SFRs are within a factor of 3 (which is 
still quite large, but our model is primitive) of the constant SFR line. We 
will define a 'star burst' to be an epoch in which the SFR has increased by a 
substantial (> 3) factor over its historical mean value. Thus short-lived 
bursts (this diagram is roost sensitive to bursts over time intervals of < 10* 
yr) are not common in field blue galaxies, and the only obvious burst 
candidate in our sample is the extragalactic HII region IIZw40. A similar 
result has been found by Thronson and Telesco (1986), although we have used a 
slightly different approach in analyzing our data. 



log L(1R) 


Figure 3. Predicted correlations for L(B) vs. L(IR) are shown with observed 
properties of spiral and irregular galaxies. 
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The L(B) vs. L(IR) diagram is not a useful diagnostic for systems such as 
M82, where Tg is high. In these cases one can use mass-to-light ratios to 
detect anomalously luminous stellar populations. This approach was pioneered 
by Telesco and Harper (1980) and by Rieke et al. (1980), who pointed out that 
very low values of stellar mass-to-light ratios are indicative of 

overpopulations of massive stars. Using our models, we would predict M*/1*(1R) 
= 2.5 in solar units for a constant SFR, and M*/l(IR) <; 0.5 for high amplitude 
SFR bursts. Observed M*/l(IR) values in the centers of galaxies suspected of 
harboring star bursts are often < 0.1, which argues very strongly for the 
existence of bursts or some other anomaly, such as an IMF that is skewed in 
favor of massive stars (see Telesco 1985). 

Another way to look at star formation histories is shown in Figure 4. 
Galaxies in region 1 have near constant SFRs in combination with either 
moderate x(tta) or dust heating by older stars to give high values of L(IR) 
relative to optical luminosities. Area 2 contains galaxies where both L(B) 
and L(Ha) are very low, and these systems are likely to be very dusty and 
effectively opaque at visible wavelengths. The physical nature of the power 
sources in these hidden galaxies cannot be determined from optical 
observations. Area 3 contains optically thin, constant SFR galaxies, and area 
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Figure 4. Irregular and spiral galaxies are plotted on a color-color diagram 
which can be used to study star forming histories. Different numbered regions 
are discussed in the text. Galaxies with high optical depths, which are 
effectively hidden at optical wavelengths, are found in the upper right 
corner, Hll-region-like objects fall to the center right, and constant SFR 
galaxies cluster near the left-center. 
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4 shows the presence of a few low optical depth burst galaxy candidates. 
Finally Zone 5 includes galaxies which have extreme excesses of L(Ha). Very 
metal-poor galaxies with low dust-to-gas ratios could populate this region of 
the diagram, but we have no outstanding candidates for such galaxies. 


5. CONCLUSIONS AND CAUTIONS 

i) Both L(IR) and L(Ha) provide useful measures of current SFRs in 
galaxies, at least to the levels of factors of 2 or 3. It is not clear, 
however, what stellar lifetime should be associated with L(IR), and there are 
indications that the average age of the stellar populations responsible for 
L(IR) varies between galaxies. 

ii) The majority of normal blue Irrs in our sample (which is not a 
complete sample; see Hunter and Gallagher 1986) have properties consistent 
with constant SFRs over the past 2-4 Gyr. Thus the IR observations give 
results that agree with the GHT optical studies (see also Thronson and Telesco 
1986), and one cannot assume that galaxies with blue colors or strong emission 
lines are in a star burst phase. Furthermore, if the SFR is roughly constant 
over time, luminosities arising from different age stellar population 
components are correlated, and thus it is not easy to isolate the stellar 
population component responsible for L(IR). 

iii) Optical depths may play a substantial role in determining properties 
of FIR-luminous galaxies. In general we expect that as the internal optical 
depth increases, L(IR) rises while L(Ha) and L(B) decline. We see some 
evidence for this trend to occur among the luminous Irrs (see Figure 4) which 
suggests that optical SFR measurements do not give a complete picture in these 
systems. 

iv) SFR bursts are not common among the dwarf and normal Irr galaxies in 
our sample. IIZw40 has all the properties that would be expected for a high 
amplitude SFR burst in a small system. The very high values of L(IR)/L(B) and 
L(IR)/L(Ha) in combination with the substantial optical luminosities of the 
luminous Irrs suggest that many of these galaxies could also be in SFR burst 
evolutionary phases. We suspect, on the basis of optical observations, that 
the majority of the luminous Irrs are results of enhanced star forming 
activity associated with galaxy-galaxy interactions. High levels of star 
formation are also seen in the strongly interacting galaxy sample being 
investigated by H. Bushouse ( 1986a, b) and interactions thus may be the main 
source of star bursts among more luminous galaxies in nearby samples. 
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DISCUSSION 

YOUNG: 

In your plot of Lh„ vs. Lb was NGC 1569 included? 

GALLAGHER: 

Yes, it looks normal. This is probably an indication that the epoch of enhanced star formation has been 
going on for ~ 10 9 yr. For time scales of 10 9 yr, the Lb/Lh q ratios rapidly approach their equilibrium 
values for constant SFR systems. 
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